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We have investigated the lowest-energy structures and elec- 
tronic properties of the Au„(n = 2 — 20) clusters based on 
density functional theory (DFT) with local density approxi- 
mation. The small Au„ clusters adopt planar structures up to 
n — 6. Tabular cage structures are preferred in the range of 
n — 10 — 14 and a structural transition from tabular cage-like 
structure to compact near-spherical structure is found around 
n — 15. The most stable configurations obtained for Aui3 and 
Auig clusters are amorphous instead of icosahedral or fcc- 
like, while the electronic density of states sensitively depend 
on the cluster geometry. Dramatic odd-even alternative be- 
haviors are obtained in the relative stability, HOMO-LUMO 
gaps and ionization potentials of gold clusters. The size evo- 
lution of electronic properties is discussed and the theoretical 
ionization potentials of Au„ clusters compare well with ex- 
periments. 
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I. INTRODUCTION 

Gold clusters are currently attracting interest as 
the building blacks of novel nanostructured materi- 
als and devices.Erla During the past two decades, the 
structures of rg<?W clusters have been— .studied both 
experimentallyLTO and theoretically.EifE 2 ] Experiments 
suggest that gold nanoclusters with diameters of 1 — 2 
nm are amorphous Im Theoretical results from empirical 
molecular dynamics simulations or ■fir.st.Tprinciples calcu- 
lations also support this argument .EEtEZroo 

In the past decade, the structure and electronic prop- 
erties of gold clusters have been intensively studied 
with various theoretical methods. Hakkinen et al. in- 
vestigated the neutral and anions of_Au2-io clusters 
with local-spin-density approximation! 2 -!! Gronbech et 
al. compared Au2 to Aus with .spin-polarized Becke- 
Lee- Yang-Parr (BLYP) functionalEj. Bravo-Perez et al. 
investigated small gold clusters upio 6 atoms at ab ini- 
tio Hartree-Fock (HF) and post-HFo level. For the large 
clusters with up to 147 atoms, Haberlen et al. performed 
a scalar relativistic all-electron density functional calcu- 
lations on several magic-sized clusters with icosahedral, 
octahedral and cuboctahedral structuresE_3. By combin- 
ing pseudopotential DFT calculations with an empiri- 
cal genetic algorithm, Garzdn found amorphous struc- 
tures in Au„ (n = 38, 55, 75) clusters.Ej In addition, 



there are some other works on the global minimal struc- 
tures of mcdjam^ized gold clusters based on empirical 
potentialsBfeJeZI 

Despite the achieved progress, there are still many 
open questions for gold clusters. For example, no direct 
experimental information is available on the structures of 
smaller Au rl clusters down to 38 atoms. Thus, accurate 
first principles calculation is important for understat- 
ing the structural and electronic properties of those clus- 
ters. Most previous ab initio calculations on gold clusters 
are limited by the presumed symmetric constraints B'cB 
A unconstrained global search on the cluster potential 
energy surface is needed.E3 In this work, we generate a 
number of structural isomers from empirical genetic algo- 
rithm simulations. These metastable isomers are further 
optimized at DFT level to determine the lowest-energy 
structures. The relative stabilities, electronic density of 
states, HOMO-LUMO gap and ionization potentials of 
gold clusters are investigated. 



II. COMPUTATIONAL METHOD 

In this work, self-consistent field (SCF) electronic 
structure calculations on gold clusters are performed by 
using a DFT-based DMol packageEa. A relativistic ef- 
fective core potential (ECP) and a double numerical 
basis including <i-polarization function are used in the 
calculations. The electron density functional is treated 
by the local density approximation (LDA) with the 
exchange-parrelation potential parameterized by Perdew 
and Wangpl. SCF calculations are carried out with a con- 
vergence criterion of 1CP 6 a.u. on the total energy and 
electron density. Geometry optimizations are performed 
with the Broyden-Fletcher-Goldfarb-Shanno (BFGS) al- 
gorithm 

the gradient and displacement, and 10 -5 
energy in the geometry optimization. 

The accuracy of the current computational scheme 
has been checked by benchmark calculations on the gold 
atom and the bulk gold solid. The ionization potential 
and electron affinity obtained for gold atom from our 
calculation are 9.78 eV and 2.51 eV respectively, which 
agpee well with experimental data: 9.22 eVEj and 2.31 
eVEil. For gold solid in fee structure, we obtained a lat- 

o 

tice parameter as 4.07 A and the cohesive energy as 4.01 
eV per atom, while the experimental lattice parameter is 



We use a convergence criterion of 10 a.u on 

a.u. on the total 
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4.08 A |-jnd experimental cohesive energy is 3.81 eV 
per atomEJ. Thus, we believe the effective core potential 
and numerical basis set used in current DFT package is 
reasonably good to describe different gold systems from 
atom to clusters and solid. 

Due to the complexity of the electron configuration of 
gold atoms (5d 10 6s 1 ), simulated annealing (SA) determi- 
nation of cluster global minimal structure at DFT level is 
rather computational expensive. Alternatively, we gen- 
erate a number of low-energy structural isomersJpij.each 
cluster size by using a genetic algorithm (GA).Ej'Ljo In 
the GA simulation, we adopt molecular dynamics with 
a properly fitted tighUbinding potential to perform local 
structural relaxtionsHJ The structures obtained from em- 
pirical GA simulations are then fully optimized by LDA 
calculation to locate the global lowest-energy configura- 
tion. The essential idea is to divide the phase space into 
a number of regions and find a locally stable isomer to 
represent each region. Our previous works show that the 
combined scheme of DFT with empirical GA is a reliable 
and efficient way for modeling the-.structural properties 
of atomic clusters up to 25 atomsEl 

III. RESULTS AND DISCUSSIONS 
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FIG. 1. The lowest-energy and metastable structures for 
small Au n clusters: n = 3 — 10. 



The obtained lowest-energy structures and some 
metastable isomers for Au ra clusters are shown in Fig.l 
(n = 3 - 10) and Fig.2 (n = 11-20). For the Au 2 
dimer, the binding energy, bond length and vibration fre- 

o 

quency are obtained as 2.43 eV, 2.55 A and 173 cm -1 , 
respectively. Our current LDA results are in satisfactory 
agreement with the experimental data (Eb = 2.30eV, 

R e = 2.47 A and lu = mcm^H'EI For Au 3 trimer, a 

o 

linear chain (3a) with bond length 2.67 A is about 0.04 eV 
lower in energy than the triangle stmctures (3b), which is 
in consistent with CASSCF studiesEj but on contrary to 
BLYP results.B Due to the Jahn- Teller instability, both 
the obtuse and acute triangle of AU3 are more stable than 
the equilateral triangle, while the energy differences be- 
tween these triangle isomers are very small, i.e., within 



0.01 eV. 

The lowest-energy structures of gold clusters with 4 — 6 
atoms are found to adopt planar forms. For Au4, a planar 
rhombus (4a) is about 0.21 eV in energy lower than the 
planar "Y-shaped" structure (4b) and 1.60 eV lower than 
the three-dimensional (3D) tetrahedron (4c). For Au5, 
the trapezoidal ( "W-shape" ) structure with Civ symme- 
try (5a) is more stable than the 3D trigonal bipyramid 
(5b) by 0.79 eV and the square pyramid by 0.94 eV. In the 
case of Aug, we obtain a planar triangle with D 3 h sym- 
metry (6a). It can also be understood in terms of "W- 
shaped" AU5 capping by an extra atom. The 3D configu- 
rations such as pentagonal pyramid (6b), octahedron and 
capped trigonal bipyramid are found as local minimum 
for Aug. The experimental photodetachment spectra of 
Au^ cluster also implies a planar hexagonal structuxa 
with low electron affinity and large HOMO-LUMO gap.O 
The planar equilibrium structures have been obtained for 
Au4_g from previous ab initio calculations^^ -EjX-, but 
can not be obtained from empirical simulationscJcB. It 
is worthy to note thai thepOjther small monovalent metal 
clusters such as Na„H3, AgEJ and Cuc3 also adopt similar 
planar configurations. 

The pentagonal bipyramid (7a) is the lowest-energy 
structure found for Au7, which is more stable than dis- 
torted capped octahedron by 0.03 eV, planar hexagonal 
structure (7b) by 0.16 eV, and capped octahedron struc- 
ture by 0.19 eV. The most stable configurations for Aus 
and Aug are largely distorted bicapped octahedron and 
bicapped pentagonal bipyramid, respectively. The onset 
of 3D lowest-energy structures starting from AU7 indicate 
a 2D^ 3D transition around the size of 7 atoms. Similar 
structural changes have been also found in alkali-metal 
clusters.^ 
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FIG. 2. The lowest-energy and metastable structures for 
medium-sized Au n clusters: n = 11 — 20. 

Some structural components related to icosahedron or 
cuboctahedron, such as pentagonal and hexagonal struc- 
tures are found in Au„ from n > 10. For example, 
both Auio and Aun can be taken as two interpenetrated 
pentagonal bipyramid, while the ground state structure 
of Aui2 is constituted by two hexagons and three pen- 
tagons. However, the lowest-energy structure for Aui3 is 
neither icosahedron (13c) nor cuboctahedron (13b), but 
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a amorphous configuration (13a in Fig. 2). This disor- 
dered amorphous configuration is about 1.44 eV in en- 
ergy lower than cuboctahedron and about 2.71 eV lower 
than icosahedron. From our calculations, the cuboctahe- 
dron is more stable than icosahedron, which has been pre- 
dicted by HaberlenJll We have also examined the bond 
lengths distributions of different isomers for Aui3. The 
disordered isomer shows broader bond length distribu- 
tions than those for the high symmetric structures. Sim- 
ilar results were found in Pti3 cluster from previous LDA 
calculation.^ 

Tabular cage structures are found for the Au^ clusters 
with n = 10 — 14. Especially, Aui4 can be taken as three 
interpenetrated pentagonal bipyramid. This behavior is 
quite different from those in the silver and copper clus- 
ters, whose ground state configurations are, usually near- 
spherical icosahedron-like structuresOOCJ. In the case 
of Auis , we find the characteristics belong to both tabu- 
lar cage and compact structure, implying a transition 
from tabular geometry to compact one. As n > 16, 
the compact near-spherical structures become dominant 
and the structures can be obtained by capping on dis- 
torted icosahedron. However, despite there are some 
icosahedral-like features in Auig , the double icosahedron 
is not found as lowest-energy structure. Similar to Auia, 
the amorphous structure (19a) is the most stable config- 
uration for Auig. The energy difference between amor- 
phous and fee structure (19b) is 1.95 eV and the differ- 
ence between the amorphous and the double icosahedron 
(19c) is 2.83 eV. Since the amorphous-structures are also 
found for larger Au„ (n = 38, 55, 75)t3, we suggest that 
the amorphous packing is the common structural feature 
of Au clusters. By using the same computational scheme, 
we have studied the structures properties of larger gold 
cluster within the size range n = 20 — 38. Amorphous 
structures are also found and the detailed results will be 
presented in a forthcoming paper. 
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FIG. 3. The average bond length of Au n versus cluster size. 

To further elucidate the growth trend of Au clusters 
and the size evolution of Au-Au interactions, we describe 
the average bond length as a function of cluster size 
in Fig. 3. Three distinct steps corresponding to planar 



structures, tabular structures and compact structures are 
found. The twice increases of the average bond length 
in Fig. 3 reflect the structural transitions at n — 7 and 
n = 15. The planar structure can be understood by 
the directionality of Au-Au bonds, implying some cova- 
lent feature in small Au clusters. As the cluster size in- 
creases, metallic bonding characters become important, 
which leads to the coexistence of planar structural and 
compact structural characters. Thus, the tabular struc- 
tures are preferred in the Auio-i4 clusters. As the cluster 
size further increases, metallic bonding eventually pre- 
vail. Therefore, compact structure appears as the domi- 
nating growth pattern of the medium-size gold clusters. 

The appearance of different structural characters in 
gold clusters can be understood by the interplay between 
5d and 6s electrons. In the small gold clusters, the 6s 
valence electrons are predominant, while the 5c? states 
are low-lying and contribute less to the electronic behav- 
ior. Thus, the small gold clusters should exhibit certain 
alkali- metal- like behaviors, e.g., planar ground state con- 
figuration, which can be described by s-orbital modified 
Hiickel modehcij As the cluster size increases, the sd hy- 
bridization increase and the contributions from d elec- 
trons become more important. Thus the clusters tend 
to adopt the more compact structures. However, as the 
short-range interaction between gold atoms is extraordi- 
nary strong, it favors a tabular structure in medium-sized 
range. The discussions below show that the 6s electrons 
still plays significant role in determining the electronic 
properties such as IPs, HOMO-LUMO gaps. 

Table I. Lowest-energy structures and electronic proper- 
ties of Au„ clusters. E\, (eV): bind energy per atom; A 
(eV): HOMO-LUMO gap; IP a (eV): theoretical adiabatic 
ionizatioiL-Dotentials; IP 6 (eV): experimental ionization 
potentialsci 



n 


Geometry 


E b 


A 


IP" 


IP 6 


2 


dimer 


1.22 


1.94 


9.81 


9.50 


3 


linear chain 


1.28 


2.70 


7.30 


7.50 


4 


rhombus 


1.74 


1.02 


8.34 


8.60 


5 


planar trapezoid (W-form) 


1.90 


1.51 


7.78 


8.0 


6 


planar triangle 


2.18 


2.06 


8.55 


8.80 


7 


pentagonal bipyramid 


2.13 


1.00 


7.20 


7.8 


8 


distorted bicapped octahedron 


2.30 


2.09 


8.19 


8.65 


9 


bicapped pentagonal bipyramid 


2.30 


0.97 


7.22 


7.15 


10 


tabular structure 


2.39 


1.03 


7.35 


8.2 


11 


tabular structure 


2.42 


0.86 


7.20 


7.28 


12 


tabular structure 


2.50 


0.82 


7.55 


8.15 


13 


tabular structure 


2.53 


0.63 


6.84 


7.70 


14 


tabular structure 


2.62 


1.58 


7.65 


8.00 


15 


tabular and compact structure 


2.60 


0.22 


7.04 


7.65 


16 


compact structure 


2.63 


0.44 


7.38 


7.80 


17 


compact structure 


2.69 


0.81 


7.33 


7.60 


18 


compact structure 


2.75 


0.92 


7.44 


7.85 


19 


compact structure 


2.77 


0.70 


6.82 


7.70 


20 


compact structure 


2.79 


0.77 


7.05 


7.82 
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Based on the lowest-energy structures, we now discuss 
the electronic properties of gold clusters. In this work, 
we have also calculated the adiabatic ionization poten- 
tials (IPs) from the total energy difference between the 
ground state neutral Au„ and the fully relaxed cationic 
Au+ clusters. The binding energy, the gap between high- 
est occupied molecular orbital (HOMO) and lowest occu- 
pied molecular orbital (LUMO), and ionization potentials 
for Au„ clusters, are presented in Table I along with ex- 
perimental IPsEfl. As shown in Table I, for small clusters, 
the binding energy increase rapidly with cluster size. As 
n > 8, the size dependent increase of binding energy be- 
come slower, corresponding to the metallic cohesion in 
the Au„ clusters. It is worthy noted that the binding 
energy up to AU20 is only 2.79 eV, which is about 70% 
of the bulk cohesive energy. 
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FIG. 4. Odd-even oscillations of cluster properties with 
cluster size: (a) Second differences of cluster energies 
A 2 E(n) = E(n - 1) + E(n + 1) - 2E(n) (eV); (b) 
HOMO-LUMO gaps (eV). 

In Fig. 4a and 4b, we plot the second differences of 
cluster total energies A2 E(n) — E(n + 1) + E(n — 
I) - 2E(n) and HOMO-LUMO gaps as a function of 
cluster size. Both the A 2 E(n) and HOMO-LUMO 
gap exhibit dramatic odd-even oscillations. The even- 
numbered Au n clusters are relatively more stable than 
the neighboring odd-sized ones and have larger HOMO- 
LUMO gap. The odd-ejpn oscillation behaviors were 
observed experimentallyQE3 and can be understood by 
electron pairing effect. Odd(even)-sized clusters have an 
odd(even) total number of s valence electrons and the 
HOMO is singly (doubly) occupied. The electron in a 
doubly occupied HOMO feels a stronger effective core po- 



tential since the electron screening is weaker for the elec- 
trons in the same orbital than for inner shell electrons. 
Therefore, the binding energy of a valence electron in a 
cluster of even size cluster is larger than that of odd one. 
In Fig.4b, the HOMO-LUMO gap for Au 2 , Au 6 and Au 8 
is particularly large (1.94 eV, 2.06 eV, 2.09 eV), which 
compare well with previous calculations by Hakkinen 
(1.96 eV, 2.05 eV and 2.04 eV).B For n = 10 - 20, sub- 
stantial HOMO-LUMO gap is found in the Auio, AU14 
and Au 18 clusters (1.03 eV, 1.58 eV and 0.92 eV). The 
high stability of Au 8 and Aui 8 can be understood by 
the effect. of s-electron shell, which is also found in silver 
clustcrsEJ. The extraordinary tabular configurations for 
the Au„ clusters with n = 10 — 14 might be attributed 
to the interplay between electronic and geometric effect. 

To explore the geometric effect on the electronic struc- 
ture, we compare the electronic density of states (DOS) 
of AU13 for the three isomers: amorphous; cuboctahe- 
dron; icosahedron in Fig. 5. The cluster electronic DOS 
shows remarkable structural sensitivity. The DOS for 
amorphous structures demonstrate more uniform distri- 
butions in comparison with the other two cases, because 
of the lower symmetry. Such structural dependence of 
electronic state can be used to identify the cluster ge- 
ometries with the aid of experimental spectroscopy mea- 
surements. 
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FIG. 5. Density of states (DOS) of the Aui3 clusters 
with different geometries: cuboctahedron, icosahedron, amor- 
phous. Gaussian broadening of 0.05eV is used. 

We further investigate the size evolution of electronic 
properties of gold clusters by examining the electronic 
density of states (DOS) of several representative clusters: 
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Au2, Aug, Auis- As shown in Fig. 6, in smallest clusters 
Aii2, the energy levels are discrete and d and sp peaks 
are clearly separated. The DOS of Aus is still sparse and 
discrete although the d and sp energy levels are gradually 
broadened. As the cluster size further increases, the d 
and sp levels broaden, shift and overlap with each other. 
Thus, continuous electronic bands are found in Auis- In 
a previous experimental study of Au n clusters up to 223 
atomsQ, the ultraviolet photoelectron spectra of smallest 
Au„ (n < 10 ~ 13) depends sensitively on the cluster 
size. The size evolution of photoelectron spectra for Au„ 
with n > 12 becomes more gradual and the spectra of 
Au8o is already quite similar to that of solid silver. 




t 

•i— ' 

O 

• i— I 

t/3 

c 

<u 

Q 



-10 -8 -6 -4 -2 

Energy (eV) 

FIG. 6. Density of states (DOS) of the Au2, Aus, and Aui8 
clusters. Gaussian broadening of 0.05eV is used. 

In cluster physics, the ionization potential is an impor- 
tant property that reflects the size-dependent evolution 
of electronic structure. For alkali clusters such as Na„, 
K„, the IPs converge to its bulk limit (work function 
of solid)-, linearly with n~ 1//3 (or l/R, R is the cluster 
radius)c£l. This behavior can be understood by a con- 
ducting spherical droplet (CSD) modeOEi In Fig. 7, the 
theoretical and experimental IPs of Au n are plotted as a 
function of n -1 / 3 and compared with the prediction from 
classical conducting sphere droplet (CSD) modeo. Our 
present results agree well with the experiments, while 
the CSD model can only qualitatively describe the size 
dependence of IPs. Similar to those in Fig. 4, dramatic 
even-odd alternative behavior is found in Fig. 7, where 
clusters with even number of s valence electrons have 
higher IPs than their immediate neighbors. In addition, 
particular higher IP values at the magic-sized clusters 



such as Au2, Au4, Au6, Aus, Au i4 and Auis is obtained. 
Some of the magic size (Au2, Au 8 , Aui4, Au«) can be 
associated to the occupation of electronic shelcS. The IP 
behavior foupd for Au„ clusters is very similar to those in 
Ag n clustersES. The significant deviation from the CSD 
model demonstrate that the Au„ clusters up to n = 20 
is still far from a piece of bulk metal. 




Cluster size n 



FIG. 7. Ionizjation potentials (IPs) of Au n . Squares: ex- 
perimental jdataH; circles: our DFT calculation; dashed line: 
CSD modeH 



IV. CONCLUSION 

In summary, the lowest-energy geometries, binding en- 
ergies, density of states, HOMO-LUMO gap and ioniza- 
tion potentials of Au„(n = 2 — 20) clusters have been 
obtained by LDA calculations combined with empirical 
genetic algorithm simulations. The main conclusions can 
be made as follows. (1) The structures of smallest gold 
clusters are planar and dominated by s electrons, sim- 
ilar to those of alkali-metal clusters and other coinage- 
metal clusters. The contributions of d electrons become 
more important and the structural transition from 2- 
D to 3-D takes place at the size of 7 atoms. (2) The 
electronic effect and geometrical effect simultaneously in- 
fluence the ground state configurations of medium-sized 
clusters, which leads to the tabular cage configurations 
with n = 10—14 and more compact structures with 
n > 15. (3) The most stable configurations of Aui3 
and Auig are neither icosahedron nor cuboctahedron, 
but amorphous structures. Remarkable difference in elec- 
tronic states are found between the structural isomers. 
(4) The odd-even alternation behaviors are found in the 
relative stabilities, density of states, HOMO-LUMO gap 
and ionization potentials of gold clusters. The even- 
numbered Au„ clusters show relatively higher stability. 
Bulk-like continuous electronic band is found in Auig 
while the electronic behavior of such clusters is still far 
from bulk metal. 
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